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Superradiance, the enhanced collective emission of light from a coherent ensemble of quantum 
systems, has been typically studied in atomic ensembles. In this work we study the enhanced 
emission of energy from coherent ensembles of harmonic oscillators. We show that it should be 
possible to observe harmonic oscillator superradiance in a variety of physical platforms such as 
waveguide arrays in integrated photonics and resonator arrays in circuit QED. We find general 
conditions specifying when emission is superradiant and subradiant and find that superradiant, 
subradiant and dark states take the form of multimode squeezed coherent states and highly entangled 
multimode Fock states. The intensity, two-mode correlations and fraction of quanta trapped in the 
system after decay are calculated for a range of initial states including multimode Fock, squeezed, 
coherent and thermal states. In order to explore these effects, the Law and Eberly protocol [C. K. 
Law and J. H. Eberly, Phys. Rev. Lett. 76, 1055 (1996)] is generalized to prepare highly entangled 
multimode Fock states in circuit QED. 
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I. INTRODUCTION 

A range of new and exciting quantum technologies have 
emerged recently that demonstrate quite detailed levels 
of coherent quantum control. In particular, exquisite 
quantum control has been achieved in highly versatile 
systems such as in circuit QED and integrated photonics. 
In integrated photonics these achievements include the 
creation of the first controUed-NOT (cNOT) quantum 
gate P] and the demonstration of quantum walks of cor- 
related photons P'. Furthermore, the Grover search and 
Deutsch-Jozsa quantum algorithms [^ and the prepara- 
tion of NOON, GHZ and W states have recently been ex- 
perimentally demonstrated in circuit QED [IHB]. Given 
these proof of principal experiments and improving fabri- 
cation techniques, next generation experiments have been 
looking towards many-body quantum systems [7] . 

However, many-body quantum systems are exception- 
ally difScult to study experimentally. Issues surround- 
ing initial state preparation, complex many-body evolu- 
tion and the impossibility of full tomography have lim- 
ited these experiments to all but the simplest of systems. 
From a theoretical perspective there are also difficulties 
in solving the complex equations decribing many-body 
system evolution. As most practical experimental sys- 
tems are few-body quantum systems, N < 20, it is not 
possible to study interesting condensed matter phenom- 
ena such as the Kondo effect and quantum phase tran- 
sitions in spin glasses. Given the technological ability, 
we seek to study new physics that can be explored in 
few-body experiments. 

One such phenomena is harmonic oscillator superra- 
diance. Superradiance is the enhanced emission from an 
ensemble of quantum systems due to field mediated build 



up of correlations within the ensemble. Initially studied 
in two level atomic ensembles by Dicke in 1954 [5], su- 
perradiance has recently experienced a revival of interest 
with its observation in quantum dots [9] and Bose Ein- 
stein condensates jT^ . However, superradiance has never 
been observed from continuous variable systems. 

In this paper, we show that recent advances in super- 
conducting quantum optics [TTJ [T^] and integrated quan- 
tum photonics [131 I14j provide an ideal environment for 
the practical exploration of harmonic oscillator superra- 
diance phenomena. This paper is organized as follows; 
in section [H] superradiance of harmonic oscillators and 
its various physical implementations are discussed. We 
introduce a collective basis, analogous to the Dicke ba- 
sis, and compare the initial emission intensities of atomic 
and oscillator ensembles in section IIIII In section IIVI a 
range of collective initial states are introduced and their 
preparation is considered in circuit QED and integrated 
photonics. The superradiance master equation is derived 
in section Ivj from which the emission intensity of an 
ensemble of oscillators is found. Decoherence free sub- 
spaces are also considered; we find a set of dark squeezed 
coherent states which do not decay and may be useful as 
a quantum memory for continuous variable quantum in- 
formation. In section [Vl] we develop a general criterion to 
classify any state as superradiant, subradiant or normal 
radiant. The emission intensity, fraction of trapped en- 
ergy and two-mode correlations are derived for a variety 
of initial states in section |VII[ Using the superradiance 
criterion, superradiant parameter regimes are found for 
these initial states. Finally, specific physical implementa- 
tions of harmonic oscillator superradiance in circuit QED 
and integrated photonics are considered in section [VIII| 
In the circuit QED implementation, the Law and Eberly 
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FIG. 1: Diagrams of experimental implementations (a) Inte- 
grated Photonics: a semicircular array of system waveguides, 
hj, collectively decay irreversibly into a semi-infinite linear 
array of bath waveguides, _Bfc. By altering the input state of 
light into the system waveguides, hj, we can control the rate 
at which light couples into the semi-infinite array, 13^- (b) 
Circuit QED: an ensemble of oscillators (squares) coupled to 
a strongly damped stripline resonator. These oscillators could 
be the transmon qubit in the harmonic regime, superconduct- 
ing stripline resonators or lumped element resonators. 



[15| protocol is generalized to prepare highly entangled 
multimode Fock states. This protocol is useful for tests 
of harmonic oscillator superradiance as it can be used to 
prepare many of the lower lying states of the collective 
basis. 



II. SUPERRADIANCE 

In this section we briefly review relevant experimental 
and theoretical work on superradiance and discuss several 
physical implementations where harmonic oscillator su- 
perradiance could be observed. Superradiance has been 
experimentally observed in a variety of systems including 
gases [TB] , ions [T7] , thin crystal slabs [TS] , quantum dots 
[9] and Bose Einstein condensates [10]. Furthermore, nu- 
merous theoretical works have considered superradiance 
in a number of physical systems such as nitrogen vacancy 
centers in diamond and circuit QED [ini [TU^ET] . 



Experiments performed to date typically observe the 
classic / ex iV^ intensity character of superradiance, how- 
ever, complicated many-body effects usually inhibit a de- 
tailed study of the dynamics of emission or the build-up 
of quantum correlations. Experimental difficulties in typ- 
ical superradiance studies arise from three sources: ini- 
tial state preparation, intrinsic losses and measurement. 
Initial state preparation is problematic as superradiant 
Dicke quantum states are highly entangled many-body 
quantum states. Most experiments cannot prepare such 
states and instead fire a pulse into the medium randomly 
exciting subsections of the ensemble [TSl HH] . Propaga- 
tion effects through the ensemble can also lead to delays 
and distortions of the superradiant emission. Secondly, 
dephasing and dissipative losses can arise due to colli- 
sions within the ensemble or via interactions with lattice 
phonons which break the coherence required for super- 
radiant decay. Finally, measurement is restricted to the 
observation of the emitted radiation rather than direct 
measurement of the ensemble dynamics, which makes a 
faithful reconstruction of the many-body dynamics diffi- 
cult. 

The vast majority of work on superradiance has fo- 
cused on ensembles of two level systems. There has also 
been some work on superradiance of multi-level atoms 
where more complex phenomena can be observed [16) . 
Such phenomena includes superradiant light beating, 
fluctuations between different light polarizations during 
superradiant emission and superradiant emission of suc- 
cessive pulses of different frequencies. Superradiance 
from an ensemble of harmonic oscillators was studied the- 
oretically in the 1970's "TI, TS most notably by Agarwal 
using phase space methods [20, 24 - 26^ . However, without 
a readily available experimental implementation, interest 
in superradiance of harmonic oscillators soon subsided 

m- 

We now consider experimental systems where har- 
monic oscillator superradiance could be observed. Firstly 
we require an oscillator ensemble. Such ensembles could 
include waveguides in integrated photonics (Fig. [l(a)[ ), 
superconducting qubits in the harmonic regime [27] . su- 
perconducting stripline [12j and lumped element res- 
onators [5S] (Fig. 1(b)), mechanical resonators, photonic 



band gap cavities, silicon toroidal microresonators and 
micro discs, among others. Secondly we require the en- 
semble to decay into a common Markovian bath. This 
could be achieved in one of two ways; either all oscillators 
are located in a small region and interact with a common 
bath (Fig. 1(a) ) or all oscillators interact with a common. 



strongly damped mode (Fig. 1(b) I. This requirement for 



a common bath is experimentally difficult and is the pri- 
mary reason why harmonic oscillator superradiance has 
not been observed to date. Thirdly, measurement of the 
intensity from the ensemble is required, either directly 
by observing the emission or indirectly by measuring the 
energy loss from the ensemble. Furthermore, desirable 
attributes of an experimental implementation would in- 
clude the ability to prepare the oscillators in a variety of 




gj > 0. Furthermore, the modes are assumed to be in- 
dependent with [a,a^ = 1, [6i,6M = 6ij and ah other 
commutators vanishing. This system can be easily real- 
ized in a variety of physical systems including the circuit 
QED system in Fig. |l(b)| [29 . The Hamiltonian ^ is 
diagonalized in Appendix A. 

The Hamiltonian can be rewritten using collective op- 
erators, 

H/h = uja^a + Lj{L + R) + 5^(6^ Cjv + aC]^), (2) 

where we have defined: the collective ladder operator 
C 
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L = 



L = l 



= Y.j=i9]bj/SN, with Qj = \jY.i=i9l; an opera- 
tor to quantify the position on the ladder R = Cj^Cn, 
an operator to label distinct ladders L = M — R and 
an operator that quantifies the numb er of quanta in the 
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hj modes M = ^27=1^% (see Fig. 2(b)). These op- 
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FIG. 2: (a) Star shaped geometry in the Hamiltonian ([l|. 
(b) Energy level diagram of the bosonic Dicke basis states, 
|d(i j), "1>|J). As all states on each ladder have the same degen- 
eracy vector, this vector is placed at the base of each ladder. 
The degeneracy vector at the bottom of each ladder is the sum 
of two unit vectors, dj^jj = Si+ej, where we assume eo is the 

zero vector. For example: jd{2,o),<&?) = (C'Jv)^C'ljO)® ^/^ 
and|d(i,2),$^)=(7;ctc.t|o^®iv^ 



interesting states and the ability to measure two-mode 
correlations within the ensemble during decay. In this 
paper, it is shown that the above requirements are satis- 
fied in an integrated photonics and circuit QED system. 



erators are analogous to Dicke's collective atomic op- 
erators [30], Ji 

where i = {+, - 



= E;=i'^]- and J2 = J^^+f^ + Jl 
,x,y,z} and ct* denote the individual 



L, J, ~ M, J4 



C 



Af 



qubit Pauli matrices, with J^ 
J_ - Cn and J+J- ~ ^ [SS]. 

We now introduce a multimode collective basis, the 
bosonic Dicke basis, which provides a convenient bosonic 
analogue to the the atomic Dicke basis. To this end we 
define a unitary transformation |66] of the N modes (6j), 
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(3) 



where, k — 1,2, ... ,N — 1, in addition to the previously 
defined Cn- As the transformation is unitary, it can 
be shown that [Cj,C^^,] = Jjjv {j,f e {l,2,...,iV} ) 
with all other commutators vanishing [23]. Noting M — 
'J2i=i ^ j^j i^ terms of these collective operators, it can 
be shown that [M , R] — [L, R] — 0. Therefore we can 
define the following orthonormal simultaneous eigenkets 
in terins of the eigenvalues R and L, 



III. BOSONIC DICKE BASIS 

In this section we introduce the bosonic Dicke basis 
and compare atomic and oscillator superradiance. The 
bosonic Dicke basis arises in the study of an ensemble of 
oscillators coupled to a central mode (Fig. 2(a) ). This is 



analogous to introducing the atomic Dicke basis in the 
study of the Tavis-Cummings model. The Hamiltonian 
for the star shaped oscillator system in Fig. 2(a) is. 



N 



N 



H/h = uja'^a + ^ ujh'jbj + ^ gj{a''bj + aSJ), (f) 
j=l 3=1 

where, all oscillators have the angular frequency, w, and 
the N oscillators couple to the central mode at the rates. 
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The eigenstate |dL,$f) is ("^^^ ^)-fold degenerate 
with respect to R, L and M; labeled by a degeneracy 
vector di — (mi, mi, ... ,711^-1)- For notational con- 
venience we write the degeneracy vector in terms of unit 
vectors, e^, i.e. d^ = X)fe=i "^fe^fc and denote the (iV— 1) 



zero vector by Gq . In addition in this new basis the ladder 
operators have the particularly simple form, 



R\ 



CN\dL,^L 



R\dL.<i>2 



R-l\ 



CU^L,^ 



= Vi?+l|dL,$ 



fl+i\ 



(7) 
(8) 



For clarity, the lower lying states of the bosonic Dicke 
basis are presented in Fig. |2(b)[ It is clear that R is the 
number of states from the bottom of each ladder and L is 
conserved on each ladder. The vertical axis corresponds 
to the total number of quanta in the state, (M), and it 
can be seen that (A/) = L + R. 

For example, when there are two oscillators the collec- 
tive operators are, 

1 




io/a.(0) 



(a) 



Ci = 


7^(.92&i-gi&2), 

t'2 


C2 = 


7^ (.91^1 +32^2), 

t'2 


M = 


b\bi + blb2 = ClCi + ClC2, 


R = 


g2 (5?^I^1 + 32^2^2 + .9132(^162 + &J^l)) , 


L = 


fi i)b\h+(i §2)^2 %f{b\b2+m), 

\ y2/ \ ^2/ ^2 



where, Q2 — \J ffi + 32- Furthermore, the lower lying 
states of the bosonic Dicke basis (|6]) for iV = 2 oscillators 
are, 




|eo,$°) = |0,0), 

|eo,$J) = ^(.9i|l,0)+52|0,l)), 

^2 



(9a) 
(9b) 



eo,$^) = 4(.9?|2>0)+32|0,2)+y2.gi52|l,l)),(9c) 



Q^ 



|ei,$?> = -^(32|l,0)-5i|0,l)), 

t'2 

|ei,$i) - -l^ {723152(12, 0)-|0, 2)) 

»2 

+ (3|-3?)|1,1)}, 
2ei,$g) = -!2(.g2|0,2)+5||2,0)-v^gi52|l,l)).(9f) 

»2 



(9d) 



(9e) 



In Fig. 2(b) these states are respectively, the bottom. 



first and second excited states on the left ladder; the bot- 
tom and first excited state on the second ladder from the 
left; and the bottom state on the right ladder. Further- 
more, for equal coupling rates, gi =32, the left ladder 
states (L ~ 0) can be expressed as. 



V ^ n=0 




R— n,n) 



(10) 



Given the above it is possible to compare the atomic 
and bosonic Dicke bases to understand the differences be- 
tween atomic and oscillator superradiance. The bosonic 
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FIG. 3: Comparison of the maximum possible initial inten- 
sity, /jv(0), from an ensemble of (a) Na atoms and (b) No 
oscillators with K quanta in each initial state. The atoms are 
prepared in \N/2,K - N/2) {N > K > 0) and the oscillators 
are prepared in |eo, "I'o'') {K > 0). Note the substantial differ- 
ence in scale in the vertical (intensity) axis as the oscillator 
emission is significantly larger for large K. 



Dicke states are simultaneous eigenstates of R and L, 
whereas the atomic Dicke states are simultaneous eigen- 
states of the collective operators, J^|^, m) — l{l + l)\l, m) 
and Jz\l,m) = 2m|/,m), where, N/2 > ? > |to| > [3D]. 
Noting M = L -f _R, we could have just as easily defined 
the bosonic Dicke basis as simultaneous eigenstates of 
L and M as the operator correspondences would imply. 
J^ ^ L,Jz ^ M. However, as we show in section [vj 
R is significantly more important in harmonic oscillator 
superradiance. 

Continuing the comparison of atoms and oscillators we 
note that the atomic Dicke basis has 2^ elements, while 
the bosonic Dicke basis has an infinite Hilbert space due 
to the infinite Hilbert space of each of the oscillators. In 
the atomic case there is a maximum excitation level for 
the atoms, \N/2,N/2) = |e)®^ and a finite number of 
quanta that can possibly be emitted from a finite ensem- 



ble of atoms, which does not exist in the bosonic case. 
This leads to large differences in the intensity of emission 
from either ensemble. In particular, the initial emission 
intensity for an ensemble of atoms prepared in the atomic 
Dicke state, \l,m) is [5], 



N 

In{Q) oc y^ (^,m|o-+(7^|?,m), 

= (^, m| J_|_ J_|Z,r7i), 
= (; + m)(/- m + 1). 



(11) 



For N atoms the intensity (11) is largest for the state 
\N/2, 0), where the initial intensity is /7v(0) ex f (f + 1). 
This is the origin of the In{0) oc iV^ hallmark of superra- 
diance initially discovered by Dicke [Hj. In contrast, the 
initial emission intensity for an ensemble of N oscillators 
with identical coupling rates, gj = g Vj, prepared in the 
bosonic Dicke state, |di,<I>|^) is [53], 

N 



= iV(di,$f|C'l,C'jv|di,<i>f) 
= NR. 



(12) 



As R is not bounded above like the atomic variables, 
N/2 > I > \m\ > 0, there is no fundamental upper limit 
to the emission intensity. 

Another example that illustrates the differences be- 
tween superradiance of atoms and oscillators is the com- 
parison of the maximum initial intensity from an ensem- 
ble of N atoms or oscillators with K quanta in the initial 
state with N > K. For oscillators with identical coupling 
rates the state with the greatest initial emission intensity 
is |eo, $^) with Jjv(O) ex NK. For N atoms with K exci- 
tations within the ensemble the greatest initial intensity 
occurs from the state \N/2,K — N/2), where from (11) 
In{0) oc NK + K - K'^. For i^ = 1 the intensities co- 
incide, however for K > I the oscillators always have 
a greater initial intensity. Interestingly, in the large N 
limit the intensities coincide. The initial intensities from 
\N/2,K - N/2) and |eo,$^^) are plotted in Fig. § for 
comparison. It is clear that the initial emission intensity 
for the system of harmonic oscillators is greater than or 
equal to the initial intensity from an equivalent ensemble 
of atoms. Furthermore the atoms have both a maximum 
possible intensity and a finite Hilbert space as there are 
only {N + 1) initial states. 



could be prepared in circuit QED and integrated pho- 
tonics. As will be shown in section |V] these collective 
squeezed coherent states arise in the study of harmonic 
oscillator superradiance from the master equation's de- 
pendence on the A^-th collective operator, Cn- These 
states may also have applications in multimode continu- 
ous variable quantum information |31) . 



A. Collective Squeezed Coherent States 

Here we introduce collective displacement and squeeze 
operators and expand the resulting squeezed coherent 
states in the bosonic Dicke basis. Collective coherent 
states can be generated via combinations of the collec- 
tive displacement operators. 



-DAr(a) = exp{aCl^ — a*Ci^) 
DkiP) = exp{/3Cl-l3*Ck) 
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9j 



X{D^,{^a),{U) 
(14) 



= Dl^ii- 



Gk+ 



■/3)nA^( 



j=i 



9j9k+i 

QkQk+l 



/3), 



where, k ^ N, and D'j{a) = exp(a6] —a*hj) is the single 
mode displacement operator. Collective coherent states 
can be expanded in the bosonic Dicke basis, for example, 

2 °° R 

|$2) = /)A.(a)|d^,$°)=e'^5]^|eo,<i>f),(15) 



R=0 



RV 



\K 



m. 



^,(;3)|eo,$o'^) = e-^ ^ ^tt'^^'^-'*?)' (16) 



L=0 



L! 



i^fe(/3)^A,(a)|eo,$[!) 



^I^ + ||3|^ 






'R\ JL\ 



(17) 



where the following identity was used, 

(7l.|di,<i>f)=V^^^^TT|di+e,.,$f+i) [k^N), (18) 

which can be easily verified using ([6]) . Displaced number 
states could also be generated with the collective dis- 
placement operators |32) . 

Experimentally, the collective coherent states that are 
the result of displacements of the multimode vacuum, 
|eo, $o) = |0)® ^, can be prepared via a product of single 
mode coherent states, |ai)i|a2)2 ■ • ■ \oLf^)-M. For example. 



IV. COLLECTIVE SQUEEZED COHERENT 

STATES AND STATE PREPARATION IN 

INTEGRATED PHOTONICS AND CIRCUIT QED 



In this section we define continuous variable states in 
the bosonic Dicke basis and discuss how these states 



|$^J = D^{a)D^{mf 



)N 



N 



e'lfa+f/3)i||^a-f/3)2(g)|f 
yN y2 yN yi ^ yn 

3=3 



a 



where, (f) is an overall phase. Therefore we require the 
ability to prepare each mode independently with a well 



defined phase and amplitude relationship between the 
modes. 

Although difficult to engineer, for the sake of complete- 
ness we now consider collective mode squeezed states. 
These can be generated via the single- and two-collective 
mode squeeze operators, 






exp 



S A2 



Uclf 



(19) 



exp (^ecA - ^C-jC/) , (20) 



where, ^ = re . With these operators it is possible 
to generate a variety of squeezed coherent states and 
squeezed number states |33| . For example the iV-th single 
collective mode squeezed vacuum is, 



I'^i) 



1 "" 



Vcosh 1 



Y^Xneneo,^r), (21) 



R=0 



where, Xr = (— 1 



)«^^(tanhr)^. 



B. State Preparation in Integrated Photonics 

We now consider how to prepare various initial states 
in the integrated photonics architecture. Firstly, the 
product coherent state \ai,a2, ■ ■ ■ , a^) can be prepared 
simply by driving the separate system waveguides with a 
specific amplitude and phase relationship. Injecting sin- 
gle mode squeezed light into each system mode is also fea- 
sible in current experiments |34| . Similarly, the product 
Fock state, |ni, n2, . . . tin), can be prepared by injecting 
Fock states from an external optical circuit. Further- 
more, a recent experiment has demonstrated that single 
mode arbitrary states of light up to the two-photon level 
can be prepared |35| . Using this protocol it is possible to 
create the unentangled state, 



N 



W 



j=l n=0 



(22) 



There is considerable difficulty in preparing the highly 
entangled multimode Fock states |d /,,<&£) and the col- 
lective squeezed coherent states. However there have 
been promising results. Recent progress in entangled 
multimode optical states have achieved the NOON state, 
(|iV,0) + \0,N))/V2, for A^ = 5 34J. Furthermore, there 
is a proposal for generating multimode entangled states 
in waveguide arrays including a four path-two photon en- 
tangled state [36] . 



driving the modes. There are also experimental meth- 
ods of preparing single resonator Fock states [37] , which 
could be used in parallel to create the product Fock 
state, \ni, n2, . . . un). Furthermore, there have been sev- 
eral proposals to generate entangled single mode and two 
mode states in circuit QED. High fidelity arbitrary single 
mode Fock states. 



E' 

n=0 



(23) 



have been experimentally demonstrated for up to A^i = 
10 using the Law and Eberly protocol [T5] in a circuit 
QED architecture with a phase qubit and a coplanar 
waveguide resonator [37] . The maximum Fock state that 
can be prepared, iVi, is limited primarily by decoherence. 
Improved qubit and resonator design and fabrication will 
allow larger Fock states to be prepared [STj [38] . Prepa- 
ration of continuous single mode states such as cat states 
and squeezed states have also been the subject of research 

[Si SO]. 

There have also been proposals for two-mode states in 
circuit QED, however, these remain difficult to engineer. 
The focus of the research to date has been on NOON 
states, due to their application in metrology [4j HI] . Re- 
cently, NOON states have been prepared experimentally 
for iV =: 1, 2, 3 at a fidehty of 0.76, 0.50, 0.33, respectively. 
NOON states with iV > 3 should be attainable with im- 
proved qubit coherence times [4l|42]. In addition, MOON 
states ^, (|M, 0) -I- |0, N))/V2, can be prepared with a 
similar protocol [^ . The Law and Eberly protocol was re- 
cently generalized [411 I43j to prepare arbitrary two mode 
states, 



w 



Ni N2 

EE 



(24) 



however, the fidelity of this state was not studied [68j . 
This protocol is particularly important for tests of har- 
monic oscillator superradiance as it can be used to pre- 
pare many of the lower lying states in the bosonic Dicke 
basis (see for example M). Research into three mode 
entanglement is still in its infancy, however, there are 
promising results such as a recent experiment demon- 
strating coherent control of photons between three res- 
onators 44 . High fidelity entanglement of iV > 2 modes 
in circuit QED is still an open problem. 



V. SUPERRADIANCE MASTER EQUATION 
AND DECOHERENCE FREE SUBSPACES 



C. State Preparation in Circuit QED 

We now consider state preparation in the circuit QED 
architecture. Similar to integrated photonics the prod- 
uct coherent state \ai, a2, ■■■ ^cun) can be prepared by 



In this section we consider damping of the oscillators 
in Fig. |2(a)| in order to provide a common bath for su- 
perradiant emission. It is assumed that each oscillator 
interacts with a zero temperature heat bath. We first 
state the master equation of the system. We next show 



that it is possible to prepare states of the oscillator en- 
semble that are decoherence free subspaces with respect 
to the central mode. Finally we derive the superradiance 
master equation and solve for the intensity of emission of 
the oscillator ensemble into the strongly damped central 
oscillator. 

Assuming that each oscillator in Fig. 



2(a) 



interacts 

with a zero temperature Markovian heat bath leads to 
non- unitary evolution of the system. The system dynam- 
ics are now described by the following master equation, 



coherent states can be generated through combinations 
of Dk{P), Sk{C) and Sk.i{£,) acting on the multimode vac- 
uum, |eo, $o), with k,l ^ N. Decoherence free subspaces 
have wide applications in quantum information |45l 146] 
and in particular the dark collective squeezed coherent 
states may have use as a quantum memory in continuous 
variable quantum information [4 7) . 



B. Superradiance Master Equation 



N 



P^-j^[H,p\ 



I?[a]p + y^ — I?[6j]p, (25) Returning to the master equation (26), we consider 



where, the Hamiltonian is defined in ([2]), the decay rates 
are k and 7j for the a and bj modes respectively and 
'D[A]p ^ 2ApA^ - A^Ap - pA^A. Furthermore, it is as- 
sumed that the losses from the bj modes are negligible 
on the timescale of interest, gj,K:^ 7^. This assumption 
requires that the oscillators in the ensemble have high 
quality factors. In the interaction picture the effective 
master equation is, 

p = ^liGNia^CN + aCJ^), p] + '^V[a]p, (26) 

where we have used the approximation, jj w 0. 



preparing the oscillator ensemble in a particular initial 
state which will then decay into the central a mode and 
rapidly leak out to be observed by a detector. The cen- 
tral mode acts as a common bath for the N surrounding 
oscillators. In this section we will show how this interac- 
tion with the common bath can enhance (superradiance) 
or reduce (subradiance) emission from the oscillator en- 
semble. 

We now assume the central oscillator is strongly 
damped, k 3> ^a^. This is equivalent to assuming that 
the a mode is always in the ground state ^48j i49j . There- 
fore we can adiabatically eliminate the a modes' degree's 
of freedom from the master equation ( 26 1 ^5] . After adi- 
abatic elimination, the reduced density matrix for the N 
modes is given by, 



A. Decoherence Free Subspaces 



ivr 



Pb 



■D[CN]pt, 



(29) 



We now show how decoherence free subspaces with re- 
spect to the central a mode can be constructed. Using 



"ij'' 



the collective commutation relations, [Cj,CL 
is easy to show that for any degeneracy vector d^ the 
states at the bottom of each ladder do not couple to the 
a mode , i.e.. 



C'«|dr,$")=0. 



(27) 



Using this fact it is clear that the state, 

\i^DFs) = \0)aY.'<L,<iL}\dL,'^l), (28) 



L,dr 



where X^l d |c(L,di,)P = li under the action of the 
Hamiltonian (pi) will not populate the a mode. As the 
central mode is not populated, the state, \'ipDFs)i will 
not decay into the bath under the evolution of the mas- 
ter equation ( 26 ) . Therefore the oscillator ensemble pre- 



pared in \ipDFs) is in a decoherence free subspace with 
respect to the central a mode ^5]. 

The collective operators provide a useful method of 
characterizing all decoherence free subspaces in the sys- 
tem. In the context of superradiance these states are also 
known as dark states, which are the special case of subra- 



diant states with no radiance. From (p7|) it is clear that 

Pk 



$fl. ) are dark coherent states. Dark collective squeezed 



where, we have also introduced the effective single oscil- 
lator decay rate, F — AQJ^/Nk. This is a useful timescale 
for the system because for a system of oscillators with 
identical coupling rates, gj — g^j, the effective single 
oscillator decay rate, F = ANg"^ /Nk, is the same as the 
decay rate of a single oscillator, Ag"^ /k. In Append ix [b] 
the density matrix elements of the master equation ( |29[ ), 
-P(i?,L,di)(t) = (dz„$f|pfa(t)|dL,$f), are found for an 
arbitrary pure initial state. 



Equation ( 29 1 is the superradiance master equation for 
harmonic oscillators [20l [24] . The superradiance master 
equation describes the evolution of a system of oscillators 
which decay into a common Markovian bath. We note 
that it is not necessary to have a strongly damped mode 
to create the common bath, as will be discussed in the 
integrated photonics system in section [VIII[ In order to 
compare systems of oscillators with independent baths to 
those with a common bath, we state the master equation 
that describes the decay of an ensemble of oscillators into 
individual zero temperature heat baths. 



F ^ 
Pi^^ T.^'^hPib] ~ b]b,pi - pib]b,), (30) 

where, it is assumed each oscillator decays at the rate, 
F/. Assuming uniform coupling rates, gj = g^j, the su- 



perradiance master equation can be expanded similarly, 



N 



Pb^ ^ Yl C^hPbb] - HbjPb - Pbblbj 



(31) 



jj=i 



The additional, i y^ j, terms in the superradiance master 
equation (31 1 can enhance or reduce the decay of en- 



semble states, leading to superradiance or subradiance, 
respectively. 

We now consider the intensity of emission of the oscil- 
lators into the common bath. The intensity of emission 
of quanta is the rate at which the oscillators dissipate 
quanta, 



lN{t)^-^^{M{t)). 



(32) 



Using the superradiance master equation ( 29 ) it can 
be shown. 



d_ 
di 
d_ 
di 



{M{t)) - Ti{Mpk{t)}^~NT{R{t)) (33) 
(i?(t)) - Ti{Rpk{t))^-NT{R{t)), (34) 



from which we find, {R{t)) = {R{Q))e-^^\ After substi- 
tuting these results into ( 32 ) the intensity can be found. 



lN{t)^NT{R[Q))e 



-NTt 



(35) 



Hence, the intensity of emission from any state is reduced 
to the problem of calculating (_R(0)). We note the inten- 
sity ( 35 ) is consistent with the phase space approach of 



Agarwal [24] after a suitable parameter substitution. 
It can also be shown from ( 29 ) that (i) is conserved. 



d_ 
di 



{L{t)) = Tr{Lpb} = 0. 



(36) 



Therefore we can associate the following physical mean- 
ings with the operators: (M) is the number of quanta in 
the system, (L) is the number of trapped (dark) quanta 
in the system and (i?) is the number of bright quanta 
that is free to radiate. It is interesting to note that the 
intensity is proportional to the initial number of bright 
quanta in the system, (i?(0)), not the initial energy in 
the system, (M(0)) = (L(0)) + (i?(0)). 



VI. CRITERION FOR SUPERRADIANCE IN 
HARMONIC OSCILLATORS 

In this section we classify the initial states in the en- 
semble that lead to, respectively, superradiant, normal or 
subradiant emission. In atomic ensembles the definition 
is clear cut: normal radiance has an intensity propor- 
tional to the number of atoms, N ^ superradiant states 
have an intensity greater than N and subradiant states 
have emission intensity less than N . Furthermore, the 



maximum intensity of emission from the superradiant 
states scales with 7V^ and the minimum intensity from 
subradiant dark states is zero. For ensembles of harmonic 
oscillators we develop a similar criterion. 

To resolve the classification of initial states for har- 



monic oscillators we expand the intensity ( 35 1 



lN{t) 






Ql 



-,-NTt 



{Rm 

N 

J2 g^gAblmm- (37) 



It is clear that the intensity depends on the initial correla- 
tions between the oscillators. Normal incoherent, uncor- 
related radiance from an ensemble of independent oscil- 
lators with independent baths have, (6|(0)&j(0)) = for 
all i y^ j [20j. On the other hand, coherent correlated ra- 
diation results from the correlations within an ensemble 
of harmonic oscillators as they interact with a common 
bath, i.e. (SJ (0)6^(0)) 7^ for some i ^ j. Therefore, 
assuming A'^ > 1, we can separate the intensity (35) into 
two parts PU) . 



lN{t)=I%{t)+lS{t), 



(38) 



where the uncorrelated and correlated parts are, respec- 
tively. 



I^iit) 



I^{t) 



NT 

A^r 



N 



-NTt 



Y^g^bimm, 



(39) 



N 



-Nrt 



^ 3,5, (&I (0)6,(0)). (40) 



«7^J 



We see that normal radiance comes from the uncorre- 
lated part and superradiance or subradiance from the 
correlated part. The correlated part is only nonzero for 
states where there exists (6](0)6j(0)) ^ for at least a 
single pair of modes, i ^ j. Therefore, we can classify 
initial states with a superradiance criterion: superradi- 
ance occurs for I^{t) > 0, normal radiance occurs when 
I^{t) = 0, and subradiance occurs when I^{t) < 0. 

Another useful quantity for classifying the radiance 
of initial states is the fraction of the total energy that 
is trapped in the system after emission. The frac- 
tion of energy that is dark (trapped) is given by i^ = 
{m)/{M{0)) = 1 - (i?(0))/(M(0)), where F G [0,1]. 
When i^ = all quanta are emitted from the ensemble, 
whereas when F = 1 no quanta are emitted. Using the 
superradiance criterion, we find that a normal state has 
the fraction, 



Fn = 1- 






■^ Normal radiance. 



and for any initial state the conditions hold, 

F < Fn <^ Superradiance, 
F > Fn <=> Subradiance. 



(41a) 

(41b) 
(41c) 



In the special case of uniform coupling rates, gj = .gVj, 
these conditions simplify considerably, 



Fk = 1 — 7T <=> Normal radiance. 



(42a) 





N " ^^' 


V -'""/ 


F < 1- 


- — <=> Superradiance, 


(42b) 


F > 1- 


- — <=> Subradiance. 


(42c) 



These conditions provide a very simple method of clas- 
sifying all initial states. It is interesting to note that for 
normal states with large N, most of the energy in the 
system is trapped. 

The fraction of dark energy monotonically decreases 
with the intensity, i.e, if two states, a and b, have the 
same energy, (M(0))q = (^(O))^ then, 



In > ^N ^ ^a < Fb 



(43) 



where I]^ is the initial intensity and Fi is the fraction 
of energy trapped of state i = a,b. Therefore, the max- 
imum radiance from an ensemble of N oscillators with 
(Af) quanta occurs when F = 0. Similarly the minimum 
intensity occurs for states with F = 1 <^ /Ar(f) = <=> 
(i?(0)) ~ 0. The fraction of dark energy is therefore a 
useful measure of the radiance from the ensemble. 



VII. SUPERRADIANCE FROM VARIOUS 
INITIAL STATES 

In this section we compare the emission dynamics for 
a range of interesting initial states. These states include 
bosonic Dicke basis states, multimode Fock states, ther- 
mal states, product single mode squeezed coherent states 
and collective mode squeezed coherent states. We deter- 
mine parameter regimes where these states are superra- 
diant according to the superradiance criterion. Several 
properties of interest are also calculated including the 
emission intensity, fraction of dark quanta and two-time 
correlation functions. 

Two-mode correlations are an important observable for 
superradiance as these correlations are the cause of the 
enhanced emission intensity in superradiance. These cor- 
relations could be used to differentiate between states 
that have the same emission intensity. Two-time cor- 
relation functions between two modes can be calculated 
using the the quantum regression theorem, 



Cij(t, 0) 



Ut 



Qn 



-{b\it)){hm 

-NTt\ 



(44) 



(l-e-^^^'')((C'^(0)S,(0)) 

(C';(0))(6,(0))j 
(Sl(0)6,(0))-(6l(0))(S,(0)). (45) 



We demonstrate the applicability of the bosonic Dicke 
states for calculating these quantities on several interest- 
ing initial states. 



Bosonic Dicke State 



We start by considering the simplest initial state in 
the bosonic Dicke basis, |di,<I>|^). As discussed previ- 
ously, this state is difficult to prepare in general, however, 
for a small number of oscillators and few total quanta 
((Af (0)) < 3) the lower lying states of the bosonic Dicke 
basis are currently within reach experimentally. From the 
superradiance master equation ( 29 1 it can be seen that 
jd^, 'I'lf) will decay down the L-th ladder and as i — > oo 



all population will reside in jd^, ^i). The intensity from 



the |dL,$f) state is. 



I^{t) = i?ivre- 



-NTt 



(46) 



an exponential decay proportional to the number of 
bright quanta in the system. 

Now that we have an expression for the intensity 
from any bosonic Dicke basis state we compare this 
to the intensity from an ensemble of atoms. We con- 
sider the emission from five atoms initially in the state 
||,|) = |e, e, e,e, e) and compare this to the emission 
from five oscillators in the equivalent state, |eo,$o)- ^"^ 
the atomic (bosonic) Dicke basis, this state is the fifth 
state above the ground state, ||,— |) = \g,9,g,9,g) 
'|eo,$o) = |0, 0,0, 0,0)) under the actio n of the coUec- 

shows the inten- 



4a 



five raising operator J+ (Cj^). Fig. 
sities for these initial states. As expected, the oscillators 
have a significantly larger intensity and decay much more 
rapidly than the atoms. This is because the oscillator 
states |eo,<I'^) with K = 2-5, have a significantly larger 
decay rate ex 5K, than the atomic state 1 1 , i^ — |) , which 
decay at the rate ex K{6 — K). The populations o f the 
atomic and oscillator states are also shown in Fig. |4(b)| 



and Fig. 4(c)[ It can be seen that the atomic states have 
a much larger spread in populations than the oscillators. 

Returning to the analysis of the emission of the bosonic 
Dicke basis state, |di,$|J), we find that the fraction of 
dark energy is F = L/{L + R). From this fraction we see 
that the higher the state is on the ladder, R, the smaller 
the amount of trapped photons. On the other hand the 
bottom of each ladder has R = and F = 1 and therefore 
all energy is trapped. Assuming uniform coupling rates, 
using the superradiance criterion (42) we find the state 
is superradiant for L < R{N — 1). 

The two-time correlation functions, Cij(i,0), are diffi- 
cult to compute in general as they are strongly dependent 
on the two modes (i,j) and the particular initial state 
|di, <!>£) considered. However, for states on the first lad- 
der, |eo, $^), the correlation function has a particularly 
simple form. 



c,j(i,0) = i?^e 

^N 



9i9j „-Nrt 



(47) 



indicating that the correlation between the oscillators ex- 
ponentially decays over time. 
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B. Multimode Fock State 

We now consider preparing the oscillators in the mul- 
timode Fock state, 



|ni,n2,... ,n]s) 




(48) 



As the oscillators are initially uncorrelated, we expect 
very different behavior from the multimode Fock state 
than from the bosonic Dicke basis state. The multimode 
Fock state is inconvenient to express in the bosonic Dicke 
basis, however it is still possible to calculate the intensity 
using (35), 



iN{t) - A^r 



S V^ 2 



g^ 



N 



3 = 1 



Ij'^J- 



(49) 



The fraction of dark quanta is. 



F = 1- 






3 "-J 



g^ 



N 2^7 = 1 "j 



= F, 



N, 



(50) 



indicating that the state has normal incoherent radiance. 
It is surprising that for uniform coupling rates {gj = g^j) 
the fraction is, F = 1 — -^, which is independent of the 
number of photons in each oscillator, Uj. Furthermore, 
for a large number of oscillators most of the excitation is 
dark as F = 1 — ^ — > 1. 

It is interesting to compare the intensity of the mul- 
timode Fock state (49) to that of the bosonic Dicke 
state ( 46 ) . Assuming we have K quanta in the initial 
state and uniform coupling rates {gj — gVj), the multi- 
mode Fock state has the intensity lN{t) = KTe^^^'^, 
whereas the bosonic Dicke basis state, |eo,<I'^), has 
lN{t) = KNre~'^^*, an A^-fold enhancement one would 
expect from superradiance. 

The two-time correlations between the modes are, 



FIG. 4: Comparison of the emission from N — 5 atoms 
prepared in the state ||,|) and A^ = 5 oscillators pre- 
pared in |eo,"l'o)- Intensity of emission: (a) the oscilla- 
tors ( 46 1 have a significantly larger emission intensity than 



the atoms ( C6 1 despite each state containing five quanta, 

(|:fl-^-l|.|) = {^o\M\^o) = 5. state populations: (b) 
The populations of the atomic Dic ke st ates during emission, 
Pk = (|,7<'-||p,||,is:-|) from ([C5|. (c) The populations 
of the bosonic Dicke states during emission Pk ~ {^(f\pb\^(f) 



from (B8l. It is assumed all atoms and oscillators couple to 



the strongly damped mode at the same rate g. 



-«j 






jj"-j, 



(51) 



where, 6i,j is the Kronecker delta. 



For i ^ j the corre- 
lation is initially zero and becomes more negative over 
time due to the decay inducing correlations between the 
modes. The initial unentangled product state becomes 
increasingly entangled as it decays. For equal coupling 
rates the intensity (49) and two mode correlations (51) 
are consistent with i24i. 



Incoherent Mixtures and Thermal states 



In this section we consider a product of single oscillator 
incoherent mixtures. These mixtures include a product 
of thermal states in each oscillator, which are particularly 
simple to prepare. Due to their classical nature, we do 
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not expect that these states will exhibit any superradiant 
phenomena. 

A product of single oscillator incoherent mixtures is 
described by the density matrix, 



N 

P(0) = (g)p,(0),with, 

oo 
n=0 



(52) 
(53) 



where, X^^o -^n = 1 ^^"^ ^'^ denote the average pho- 
ton number in each oscillator by fij = X]n=o"'^n- ^'-"" 



example a thermal state in the j'-th oscillator has, 



pj 



i^jY 



{fij + i)"+i ' 



(54) 



where, fij — l/(e''"' '^-^j _i) is the average thermal photon 
number of the j'-th oscillator at temperature, Tj. Fur- 
thermore the multimode Fock state considered in ( 48 1 is 
also an incoherent mixture with P^ = 6n n ■ 



As the density matrix (52 1 factorizes it is clear that 



(6{(0)6j(0)) = OVi 7^ j. This imphes that the correlated 



part of the intensity (401 is zero, indicating that an in- 
coherent mixture has normal radiance. The fraction of 
dark quanta, 



F= 1- 



V^^' 2- 



Ql 






= F, 



N, 



(55) 



supports this conclusion. Therefore, the emission from 
a system of oscillators in a thermal state is normal. Al- 
though the oscillators are initially uncorrelated, the emis- 
sion into the common bath induces correlations as can be 
seen from the two-time correlation function, 



^^,3 



(t,0) = 



9igj 
'5% 



(1 



^-Nrt 



hj+kjri]- 



(56) 



Finally, the intensity of emission from an incoherent 
mixture is similar to the multimode Fock state. 



-Nrt 



N 



iN{t) = A^r- 



Gl 



E^: 



^"^j- 



(57) 



For uniform coupling rates, gj — g^j, and average pho- 
ton numbers, fij — nyj, the intensity simplifies to. 



I^{t) ^fiNTe 



^NTt 



(58) 



We see that the incoherent mixture has the intensity one 
would expect from normal radiance, I^{t) oc N . 



of particular interest because of the close connection be- 
tween coherent states and classical harmonic oscillators. 
The emission from an ensemble of classical oscillators 
with random phase scales with /^r ex N , whereas for clas- 
sical oscillators that are in-phase the emission scales as, 
In oc A^^. The in-phase intensity is clearly reminiscent 
of superradiance. 

In addition to coherent states, we also consider the 
possibility of single mode squeezing [M] in each oscilla- 
tor, as each mode can be prepared independently. There- 
fore, we consider the initial state where the oscillators are 
prepared in a product of single mode squeezed coherent 
states. 



|ai,Ci)i|a2,C2)2---|aAr,C 



NIN 



N 



■Qga,ft]-a*6,gi(e*6^^-C,(h])^)|Q^C 



)N 



(59) 



where, ^j = tje^^ . The intensity of emission for this 
state is. 



lN{t) = 




N 



E 5i sinh^ rj 



(60) 

The most important feature of the intensity is that the 
amplitudes, aj, add coherently. This feature is most ap- 
parent for the special case of uniform coupling rates, dis- 
placements and squeeze amplitudes {gj — g, aj — a, tj — 

ryj), 

I^{t) ^T{N^\a\^ + N sinh^ r)e-^^K (61) 



We see that the intensity has the /^(i) oc A^^ charac- 
teristic of atomic superradiance. However, the intensity 
from the squeezing contribution only scales linearly with 
the number of oscillators, a characteristic of normal radi- 
ance. In fact for squeezed vacuum states (a = 0) the in- 
tensity is I%{t) ex N, whereas for coherent states (r = 0), 
I%{t) oc N^ 



In order to classify the radiance of the state (59), we 
find the fraction of dark quanta. 



^N 



N 



F = 1 



Ej=i ffjajh + Ej=i 9j sinh 



g^ 






sinh 



(62) 



Using this fraction and the superradiance criterion we 
find that the state is superradiant if ^70], 



N 

IE 



N 



(63) 



D. Squeezed Coherent States 

Another state which is straightforward to prepare is 
a product of single mode coherent states. This state is 



Therefore, the superradiance of the state is independent 
of the squeezing amplitudes, rj. If each oscillator ampli- 



N 



^N 



tude has a random phase, | J2j=i 9j^j\ — X^j^i 9jWj\ ^ 
and it can be shown that the state has normal radiance. 
This phenomenon is analogous to classical oscillators. 
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FIG. 5: Fraction of trapped quanta from A'^ = 10 oscillators 

„ie\(g) 10 



prep ared in a product squeezed coherent state |a,re' )^ 
( |64[ |. It is assumed all coupling rates, displacements and 
squeeze amplitudes are identical, Qj = g,aj — a.r-j = rVj. 
Normal radiance occurs for F^ = Fjv = 1 — ^ = 0.9, super- 
radiance occurs for F^ < 0.9. We see that large squeezing 
amplitudes trap dark quanta in the system and the maximum 
radiance occurs for coherent states (r = 0). 



In the special case of uniform coupling rates, displace- 
ments and squeeze amplitudes {gj = g,OLj = a,rj — 
ryj), the fraction of trapped energy becomes, 



F^ 



1 



a- 



N 



sinh r 



sinh^ ; 



(64) 



which is shown in Fig.pl It is clear that _F^ e [0, 1 — ■^], 

with F*^ = for r = and F'^ = 1 - ^ = i^AT when a = 
0. Therefore for |a| > the state is always superradiant, 
reaching the maximum intensity of emission when there is 
no squeezing, r — 0. Furthermore, single mode squeezing 
has the effect of trapping photons in the system as r > 
=4> F^ > 0. Finally, the squeezed vacuum state, a = 0, 
exhibits normal radiance. 

The two-time correlation function between the modes 



^«j 



(t,0) = 



9i9j (. 



^-Nrt 



) sinh^ 



-Si_j sinh rj. 



(65) 



This coin cide s with the incoherent mixture correlation 
function (56) after the replacement, sinh r,- — >■ n,. We 



note that for no squeezing, tj = 0, the correlation func- 
tion is zero, Cij{t, 0) = 0. Therefore correlations between 
the modes are not induced during the decay of product 
coherent states. This is surprising because in atomic su- 
perradiance, superradiant states are highly entangled and 
there is a large correlation between emitters during emis- 
sion. 

Summarizing the above results, we see that a prod- 
uct of single mode coherent states can be superradiant 
(63), whereas a product of single mode squeezed states 



possesses many of the properties of an incoherent mix- 
ture. For example, consider the state with no displace- 
ment (aj = 0), |Ci)i|'C2)2 ■ • • |CAf)w- The average photon 



number in each mode, 
the replacement sinh ? 



'J' 
-^ n 



is {b]b,) 



sinh 



rj . Under 



J' 



we find that the inten- 



sity ( 60 ) , fraction of trapped quanta ( 62 ) and two-time 



correlation function ( 65 ) all reduce to the corresponding 



equations for the incoherent mixture, (57), (55) and (56) 
respectively. Furthermore from the fraction of trapped 
quanta we see that the radiance is normal. We conclude 
that single mode squeezing does not lead to superradiant 
behavior. 



E. Collective Squeezed Coherent States 

In this section we consider preparing the oscillator en- 
semble in the collective squeezed coherent states from 
section |IV[ The intensity, fraction of dark quanta and 
conditions for superradiance are listed in Tablelllfor three 
particular initial states: a coherent state in the iV-th col- 
lective mode, |<&2)' ^ coherent state in the fc-th collec- 
tive mode, 1$? ), (fc 7^ N) and a product of collective 
squeezed coherent states. With reference to Fig. |2(b)[ 



I $2) ^^'^ coherent states in the bosonic Dicke basis that 



occupy states vertically on a single ladder denoted by, 
di. On the other hand the coherent states, |<&n',), pop- 
ulate bosonic Dickc basis states along a diagonal: start- 
ing from the 'ground' state 180,$^), the state occupies 
one position on the i?-th rung on each ladder labeled by 
di = Lek- Furthermore the product collective squeezed 
coherent states in Table 11] with aj ^ OVj will have some 
population on all rungs of all ladders. 

We first discuss the collective coherent state, |<&2)- ^^ 
the decay does not affect the ladder label, evolution un- 
der the superradiance master equation ( 29 ) will drive the 
population down the ladder to the dark state |di,$°). 
Therefore the intensity is independent of the undisplaced 
initial state |di,<I>2^) [7T]. Superradiance only occurs if 
|ap is large with respect to the number of dark quanta, 
L. Two-time correlation functions between the modes are 
in general difficult to calculate as it depends on the par- 
ticular degeneracy vector, d/,. However in the simplest 
case. 



N 



\n) = 



9j_ 

Qn 



131 



(66) 



the state is a sp ecial case of the product coherent state 

9i 



m section 



VII D 



with (Xj — ct-d^ and S^j — 0. From (65) 



we find that no correlations are induced during emission 



as c. 



*j 



(t,0)=0. 



Next we consider the collective coherent state 1$? ). 

As the displacement operator Dfc(/3) does not add any 
bright quanta to the system the intensity coincides with 
that of the bosonic Dicke basis state (46). This state 
is only superradiant for small displacements, /3. Similar 
to the state |$2) '^^ note that the two-time correlation 
functions depend on the specific degeneracy vector, d^, 
however, for the simplest case, | <&«,), it can be shown 
that Cjj(i,0) =0. 
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TABLE I: Collective Squeezed Coherent States 



State Intensity (/jv(i)) 

nr.ii'.(«i)4(G)|eo,-l'EJ> TVrdaivP +sinhV^)e-^" 



Condition for Superradiance 
|Qp(iV- 1) > L " 
i?(iV-l) > i/?p " 

I |2 I • 1 2 ^ 1 v^W 2- 

lajvl +smh tn > g^ l^j^^gjU 



Fraction of Dark Quanta [F) 

L 

(L+|q|2) 

1 l^ivl +sinh r^v" 



"Assuming uniform coupling rates, gj = g Vj 

''(k^ N) 

'^fij = |ajp + sinh^ Tj 



Finally, we consider the product collective coherent 
state in Table [l] As only the A^-th collective mode de- 
cays, displacement and squeezing of the other modes has 
no effect on the intensity. It is interesting to note that for 



uniform displacements and squeezing, aj — a, 
the state has normal radiance, F — 1 



= F^ 



rVj, 
As 



discussed above the two-time correlation function is dif- 
fucult to compute, however for the A^-th collective mode 



squeezed vacuum (21), \^q), it can be found, 



'^'■•j 



(i,0) = |f^e-^i^*sinh2 



rN- 



(67) 



interaction between the modes in the ensemble during 
superradiance. Firstly, low fidelity state preparation can 
lead to transitions along undesired emission pathways, 
which can reduce the intensity of emission by populating 
dark states. Secondly, losses such as those from low-Q 
resonators are a form of non-radiative decay and also lead 
to undesired transition pathways, which in turn results 
in reduced intensity. Lastly, low efficiency detectors and 
dark counts can also reduce and/or distort the measured 
emission intensity. The integrated photonics and circuit 
QED systems we consider satisfy the above requirements. 



^N 



The A'"-th collective mode squeezed vacuum, \^q), pos- 
sesses many of the properties of the superradiant bosonic 
basis state, |eo,$o^); from Table K^ the intensity, fraction 
of trapped quanta and two-time correl atio n functio n of 
|$q), are identical to those of |eo 



$^) 



(|46|), F^O, (|47D 
R. 



respectively, under the replacement sinh r^ — > R. This 
is in contrast to the single mode squeezed states in section 
11 D which were similar to normal radiant incoherent 
mixtures. We conclude that squeezing of the A^-th col- 
lective mode can lead to superradiant behavior. 



VIII. IMPLEMENTATIONS 

In this section we examine two experimental imple- 
mentations which could be used to observe superradiance 
of harmonic oscillators from a variety of initial states. 
We consider both optical and microwave implementa- 
tions. In the optical regime, we discuss direct write inte- 
grated photonic waveguide arrays in integrated photonics 
[Ol UM [501 [5T] . In the microwave regime setups based 
on the circuit QED architecture are considered including 
stripline [52j and lumped element resonators [28j and the 
transmon superconducting qubit tuned to the harmonic 
regime [27] . In each system initial state preparation is 
considered and it is shown how the superradiance mas- 



A. Integrated Photonics 

Here we consider an integrated photonics implementa- 
tion as shown in Fig. 6(a) The implementation consists 



of two sections; one for state preparation and the other 
for the observation of superradiance. The Hamiltonian 
for the state preparation section is, 

N N-1 

H/h = ojY, % + J ^ {byb,+, + b]^J>,), (68) 

where, J is the coupling rate between system waveguides. 
For an initially unentangled input state, evolution under 
this Hamiltonian can entangle the modes. For example, 
for the input state |V'(0)) = |^i)^2, • ■ • ^n^) the evolved 
state is [S3], 



^ 1 

w)) = n ^ 



J = l 




|0)^ 



)N 



(69) 



where 



ter equation (29) describes the dynamics under a specific ■^j,q — 



parameter regime. 

In choosing an experimental implementation there are 
three main requirements: high fidelity state preparation 
of a variety of states, small losses and high fidelity mea- 
surement of the superradiant emission. Two mode cor- 
relation measurements are also desirable to measure the 



A^- 



N 



-i(c^-(-2Jcos(J^)t)) 



sin( 



qpn 

N + 1 



sm( ). 

^A^ + r 



We now show how superradiant, normal and subradi- 
ant states can be generated in a three waveguide system 
by injecting a single photon into one of the modes. Us- 
ing (69) we find that the states, |V'i(0)) = |1,0,0) and 
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IIW 



-B4 -B5 



results in a subradiant state and injecting a single photon 
into the second waveguide results in a superradiant state. 
Normal states can be prepared in either case for t — 
mT/V2J. At t ^ t* := {n + ^)n/V2J, the first state 
is most subradiant {Fi = 5/6) and the second state is 
most superradiant (F2 = 1/3). Therefore we suggest 
engineering the length of the waveguides such that the 
Hamiltonian ( 68 ) evolves for t = t* . We note that more 



(a) 



Strongly 
f Damped 



Superradiant 
Emission 



Preparation / 




Tunable 
Couplers 



Tunable 
Couplers 

flj(*) 



(b) 



FIG. 6: Implementations with initial state preparation a) 
Integrated Photonics b) Circuit QED. Integrated Photonics: 
Three system modes (left) initially coupled to one another fan 
out and collectively couple to a linear array of bath waveg- 
uides (right). Circuit QED: strongly damped a mode for ob- 
servation of superradiant emission (top) connected to tunable 
couplers (squares) which in turn connect to the system modes 
(circles). The system modes are coupled to a qubit via an- 
other set of tunable couplers for initial state preparation. 



|V'2(0)) = |0, 1, 0), wiU evolve into, 

\Mt)) - a;((z + 1)2|1, 0,0)- 72(^2 „ 1)10,1,0) 

+ (z-l)2|0,0,l)), (70) 

\Mt)) = e-*-*(--^sin(y2Ji)(|l,0,0)-f|0,0,l)) 



+ cos(^Ji)|0,l,0)), 



(71) 



where, x — ie^*("+v2./)t ^nd z — e*^^"^*. Assuming 
uniform coupling rates when these states enter the su- 
perradiant section of the device {gj ~ g^j) the fraction 
of dark energy for these states are. 



Fi = 



F2 = 



1 

12 
1 

6 



cos(2\/2Jt) > Fat, (72) 

cos(2V2Ji) < Fn. (73) 



where, F/v = |. For t ^ rnr/y^J, (n S Z+) it is clear 
that injecting a single photon into the first waveguide 



complex states can be engineered by injecting more than 
one photon. 

After the preparation section, the photons enter the 
second section for the observation of superradiance. The 
superradiance section consists of a semicircular array of 
system waveguides, bj, surrounding the end of a semi- 



infinite equispaced array of bath waveguides, Bk- 
Hamiltonian for this system is. 



The 



N 



Oscillator 
Ensemble 

b, 



H/h ^ LoY, f>ybj + GNi&NBi + B\Cn) 

+ uY. ^l^'^ + f T.^bIb,+^ + Bl^.B,) (74) 



fe=i 



fe=i 



where, each waveguide has the same angular frequency 
uj, the system waveguides couple to the first bath waveg- 
uide, Bi, at the rates, gj, and bath waveguides nearest- 
neighbor couple at the rate, k/2. By design, dissipation 
within the waveguides is minimal and can be ignored on 
the timescale of interest [131 ISOl IM] • 

We make two fundamental assumptions; 1) that the 
system waveguides do not couple to one another and 
2) that the bath waveguides couple to one another at 
a larger rate than the system modes couple to the bath, 
K ^ Gn- Physically the case of TV = 3 syst em waveg- 
uides is relatively easy to envisage (Fig. 6(a)), however, 
additional system waveguides are problematic [72 . The 
second assumption requires that the distance between the 
system waveguides and the first bath waveguide is large 
compared to the distance between each bath waveguide. 
Starting from Hamiltonian ( 74 ) , we can trace over the 



bath degrees of freedom to obtain the effective master 
equation of the system [JS]. It can be shown that the 
effective master equation of the system is the superradi- 
ance master equation ( 29 1 [55l [S^ . Experimentally the 
emission intensity, /^v (t), is found as the sum of the in- 
tensity in each bath waveguide. Furthermore, correlation 
functions can also be measured in this system [2]. 



B. Circuit QED 

We now consider harmonic oscillator superradiance in 
a circuit QED system. We first introduce the system, 
then show how the dynamics are well described by the 
superradiance master equation and finally describe how 
various initial states can be prepared in the system. The 
circuit QED system, depicted in Fig.|6(b)] consists of two 
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sections, the top section is for the observation of super- 
radiance and the bottom section is for state preparation. 
The system is highly versatile as many different kinds 
of oscillators can be used, including stripline [52] and 
lumped element resonators ,28] and the transmon qubit 
operated in the harmonic parameter regime [27] ■ The 
unitary dynamics of the circuit QED system is described 
by the Hamiltonian, 



N 



"t"^ Y-?t£ ^^W 



H/h = uja'^a + ujy^^b\b^ 

N 



&,+£{t)a++£*{t)a^ 



N 



J2~9nt)ib]a + a%) + ^gj-(i)(S]<7_ + S,<7+),(75) 



where from Fig. |6(b)[ the top oscillator represents the 
a mode, which couples to the N central oscillators, 6j, 
via the upper set of tunable couplers at the rates g'j{t) 

[STlllH]. The central oscillators, 6j, also couple to a su- 
perconducting qubit via lower set tunable couplers at the 
rates, g\{t). The qubit is driven at the complex ampli- 
tude, £{t), with tunable transition frequency, 17 (i). Dis- 
sipative dynamics are described by the master equation, 



N 






V[b,\p 



(76) 



where, k is the decay rate of the strongly damped a mode 
and the collective modes hj decay at the rates, 7j. Fur- 
thermore, the qubit undergoes energy relaxation at the 
rate 7| and dephases at 7^. We assume that the dissi- 
pative rates, 7j,7|, and 7^, are small compared to all 
other system parameters. This assumption is required 
to achieve high fidelity state preparation and in order to 
approximate the system dynamics by the superradiance 
master equation. 

To observe superradiance it is necessary to prepare an 
initial state in the oscillators as discussed below, then de- 
couple the qubit by either tuning the qubit 's transition 
frequency far from a; or setting the lower tunable cou- 
plers to the off position, g''At) = 0. The system dynamics 



are therefore described by the master equation in ( 26 1 . 
Assuming the d oscillator is strongly damped, k 3> QN^ 
we can adiabatically eliminate the a degrees of freedom 
as in section JV] and arrive at the desired superradiance 
master equation ( 29 ) . Therefore harmonic oscillator su- 



perradiance can be observed by preparing the oscillator 
ensemble in a superradiant state and measuring the emis- 
sion of photons into the a mode. 

We now describe a method to prepare highly entan- 
gled multimode states using the system in Fig. |6(b)[ For 
preparation, the a oscillator is decoupled by tuning the 
top row of couplers to their off state, g"(i) = 0. We 
consider a generalization of the Law and Eberly protocol 
to prepare entangled states of several modes |T5]. The 



Law and Eberly protocol prepares an arbitrary single 
mode state (231 by successively exciting the qubit into 



a superposition of the ground and excited states then 
partially transferring quanta from the qubit to the os- 
cillator. The protocol requires the ability to switch be- 
tween two interactions: 1) classical driving of the qubit 
for single qubit rotations and 2) a Jaynes Cummings in- 
teraction between the qubit and the mode. Furthermore, 
brief pauses are required between these two interactions 
to adjust the phase of the excited state of the qubit rel- 
ative to the ground state |37| ST] . By varying either the 
coupling rates or the length of interaction time intervals, 
Law and Eberly demonstrated that an arbitrary single 
mode state ( 23 ) could be prepared [H] . 



We now discuss how each interaction can be generated 
in the circuit QED system in Fig. |6(b)[ The first inter- 
action we require is the classical drive for qubit rotation. 
In the interaction picture, qubit rotations are generated 
by evolution under the Hamiltonian, 



Hd = Eg4 



E* 



(77) 



where the drive is set to a constant amplitude, £{t) = 
£. This Hamiltonian can be obtained from (75) by two 
different methods. Either the tunable couplers could be 
switched off (g' (t) = 0) or the qubit transition frequency 
could be tuned far out of resonance, Q,{t) ^ w or il{t) <C 
w [37]. 

The second interaction we require is the Jaynes Cum- 
mings interaction. In order to attain the Jaynes Cum- 
mings interaction, the qubit is tuned into resonance with 
the modes, i.e fl(t) = w, and the tunable couplers are 
tuned to a predefined constant rate gUt) = gj- In the 
interaction picture the Hamiltonian ( |75[ ) now describes 
the Jaynes Cummings interaction, 



Hjc = g{d^^- +do-+), 



(78) 



where, d = J2j=i9jbj/g, g = yJ2j=i9'j and [d,d''] = 1. 

We note that d can represent a variety of modes, includ- 
ing Cj, where j = 1, 2, . . . , TV. The use of the collective 
mode, d, rather than a single mode, bj, is key to our 
proposal. 

Using the Law and Eberly protocol, sequentially alter- 
nating between evolutions under Hd and Hjc, with brief 
pauses to correct the phase between the ground and ex- 
cited state of the qubit, it is possible to prepare the state, 



TVi 



{dry 



|V^) = 5:c„^|0)«^ 



(79) 



are arbitrary complex amplitudes with 
= 1 and iVi is the maximum Fock state, 
which is limited due to dissipation. For example, when 

d — Cn, the state, J2rLo'^r\^0'^o)j '^^^ t)e prepared. 
For large iVi and small |^| this state could approximate 

the collective squeezed vacuum |$o)- Similarly, when d = 



where. 
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Ck {k ^ N), we can prepare the state, X^lIo '^^1-^^*;' *^fc)- 
Furthermore, for d = {Ck + Cn)/V^, {k ^ N) and 



= 5n 2 the state, 



V2 



\0)' 



)N 



2^2 



V2|eo, $^) + \/2|2efc, $«> + 2|efc, $} 



can be prepared. 

We have shown that a variety of muhimode entangled 
states can be prepared using the circuit QED system in 
Fig. |6(b)[ These states are useful in the study of har- 
monic oscillator superradiance as they include several of 
the collective states studied in section |VIIJ In addition, 
multimode states prepared by the generalized Law and 
Eberly protocol could be useful for metrology [55] . 



IX. CONCLUSION 

In summary, we have studied superradiance in an en- 
semble of harmonic oscillators. We have shown that har- 
monic oscillator superradiance should be observable in 
a variety of systems including waveguide arrays and su- 
perconducting resonators. The bosonic Dicke basis was 
introduced which proved very convenient in calculations, 
particularly in classifying dark states and in analyzing 
the dynamics of the system. We also defined a superra- 
diance criterion, which can be used to classify arbitrary 
harmonic oscillator ensemble states as superradiant, nor- 
mal or subradiant. 

Furthermore we characterized a range of states includ- 
ing bosonic Dicke basis states, multimode Fock states, in- 
coherent mixtures, thermal states, single mode squeezed 
coherent states, collective coherent states and collective 
squeezed coherent states. For each state the emission in- 
tensity, fraction of dark energy and two-time correlation 
functions were calculated. One of the most interesting 
findings was that single mode squeezing has normal radi- 
ance, whereas squeezing in the A'^-th collective mode was 
superradiant. Finally we generalized the Law and Eberly 
protocol to create highly entangled multimode states in 
circuit QED. The protocol should allow the preparation 
of many of the states discussed. 

We anticipate our results to be useful in further stud- 
ies of harmonic oscillator superradiance and other multi- 
mode phenomena. The systems we considered may have 
applications in continuous variable quantum computing 
and the decoherence free subspaces may be useful as 
quantum memories. In addition, the state preparation 
protocol may produce states useful for metrology. Fu- 
ture work will consider entanglement and the effect of 
dissipation on the dynamics in this system. 
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Appendix A: Eigenstates 

In this section we diagonalize the Hamiltonian (pi) in 
the bosonic Dicke basis introduced in section llTTl In order 
to find the eigenstates we split the Hamiltonian into two 
commuting parts. 



where, 

H' = uja^'a + ujR + QN{a^CN + aCj^). 
We now introduce the normal modes, 
1 



d± 



V2 



(a±Cjv) 



(Al) 
(A2) 

(A3) 



where, [di^d,] 
is diagonal. 



(5ij. In terms of these normal modes H' 



H' ^{u + gN)dld+ + (a; - gN)dld^ (A4) 



Therefore the eigenstates of the Hamiltonian (Al) are 
simultaneous eigenkets of L and H', 



H\l, n+, n_ 



£^« + ,n_l^"+'^i- 



(A5) 



where, £' „_ = uj{l + n+ + ?i-) + GN{n+ — n_) and 
n± = {dj-d±). The eigenstates are explicitly. 



^,„^,„_^^(Ap^,„.^^, 



•J' 



n+\n 



Ir7 I 



(A6) 



Appendix B: Density Matrix Elements of the 

Superradiance Master Equation for Pure Initial 

States 

We now solve for the matrix elements and emission 



intensity of ( 29 ) given an arbitrary pure initial state of 



the bi modes, 



\i^m 



L,R=0 di. 



a(fl.L,d,)|dL,'i>f), (Bl) 



where, J2ZR=oT.di^\"'iR-L,dL)\'^ = 1- To proceed 
we introduce the matrix elements, P(r,l.Al){^) — 
{dL,<^1\p}j[t)\dL, $£). In terms of these matrix elements 
the superradiance master equation is transformed into 
the difference differential equation. 



d 
dt 



^,P{R,L^A,){t) ^NT{{R+ l)P(«+i,L,d,) W - RPiR.L.d,){t)) 

(B2) 
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We can re-express the difference equation (B2| using 
vector notation, 



P(t) = NTAP{t), 



(B3) 



where,P(i) = (P(o,L,d,) W, P(i,L,d,)(t),P(2.L,d,) W, • ■ • )^, 
P„(i) == P{n-i,L,dL){t), and, 

(1 - i) if i = j 
otherwise. 



Appendix C: Superradiance of Atoms 

In order to compare the superradiance of oscillators 
to the superradiance of atoms, we briefly review some 
relevant results [El [61]. Consider a system of N two- 
level atoms coupled to a strongly damped mode, a. In the 
interaction picture the system evolves under the Tavis- 
Cummings Hamiltonian, 



H = g{J-a^ + aJ+), 



(CI) 



To solve the system of ODE's ( B3 1 we diagonalize A using 
the Pascal matrix [60 , 



otherwise. 



After some manipulation it can be shown that, A — 
[B-^YDB'^ where, Aj = (5*^(1 - i)- Therefore the 
eigenvalues of A are A; = (1 — i) and the eigenvectors, 
Vi^ are the columns of (B^^)'^ . Hence for a given vector 
of initial conditions P° = (|Q(ox.dz.)lMQ(i.L dz,)P, ■ ■ ■ )^, 
where P^ = |a(fc-i,L,di,)P i the solution to (B3) is. 



m - E 



Ce^'^^S, 



(B4) 



4=1 



where, the coefficients c = (ci, C2, . . . )'^ can be deter- 
mined from the initial condition, 

c = B^P°. (B5) 

We can rewrite the n-th element of this solution as, 

'i~ 1\ /k- V 

i—n k—1 
oo 
y^ pOg(l-n)A^ri/^ _ ^-Nrt\k-n i 



oo CX3 



P"W = EEp"^'"''"(-ir 



n — 1 J \i — 1 
k-1 



k—n 



n — 1 



(B6) 



where, P„(t) = -P(„_i,L,dz.)(0 and P^ ^ (|a(fc-i,L,d,,)P- 
This result can also be written in terms of the matrix 
elements (instead of using vector notation). 



Pi 



)W = E 



1 2 -B.NTt 

(fl,L,di)ltj = 2^ \a(k,L,dL)\ e 

k=R 



Nrt\k-R 



il-e-"'') 



R 



(B7) 



For example for the single bosonic Dicke basis state, 

|V'(o)) = |d,.,$f), 



iR,L,dL 



.<')^a)' 



'(l-e-^"^i 



(B8) 



The populations for five oscillators prepared in |eo,$o) 



are displayed in Fig. 4(c) 



where, all atoms couple to the field at the rate, g, and 
the collective operators are J± — X]i=i ^i ■ The master 
equation for this system is, 



P = -i[H, p] + '^v[a]p + E y ^['^7]/^ + y^t'^J] 



where 7^ and 7^ are for the jth atom the spontaneous 
emission and dephasing rate respectively and k, is the 
field decay rate. Similar to section |V) we assume n ^ 
g ^ li'l^ and adiabatically eliminate the field degrees 
of freedom to arrive at the superradiance master equation 

mm, 






(C2) 



Expanding the superradiance master equation in the 
atomic Dicke basis, we obtain the following equation for 
the matrix elements, 



P{l,m,t) = T[{l-m){l + m+l)P{l,m + l,t) 

- {l + m){l-m + l)P{l,m,t)]. (C3) 

where, P{l,m,t) = {l,m\pq{t)\l,m) and F = -^ is the 
single atom decay rate. Assuming the atoms arc initially 
aU excited, P{l,m,0) = P{N/2, N/2,0), and noting I is 
conserved by the superradiance master equation ( C2 ) , we 
can make a change of variable to n = / — tti = 0, 1, ..., iV. 
Here, n, corresponds to the number of photons emitted 



from the ensemble. In terms of n (C3) becomes 



P{n,T) = {N -n + l)nP{n-l,T)-{N -n){n + l)P{n,T), 

where, we have rescaled the time t ~ Ft. Using the 
Laplace transform the matrix elements are found to be. 



1 -rV {N-i + l)l 

-s + N ^i- s + (N -~i)(i + l)- 



Pi-^-) = ^-'{vir^I[^^'i^r^^^^^^.}- (C4) 
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where, C ""^{/(s)} — f{t), is the inverse Laplace trans- 
form. For example for A^ = 5 the populations are, 



P(0,r) 



-5r 



be- 



-5t 



3 3 

4«^"" -lOe- 



lOe- 



(C5a) 
(C5b) 

(C5c) 



- UOe'^W (C5d) 



P(4,t) = - 



220e- 



ISOe" 



320e- 



3 3 

320e-»^T + 80e~^^r (C5e) 

9r 125e-8- 



P(5,t) = l-100e^^^ + 



3 



172e-5^ 



These populations are plotted in Fig. |4(b)[ where we 
use the substitution, n = 5 — K. The emission inten- 
sity can be found from the populations via /jvIt) = 

57 X]n=o "-P("7 ''')• For example for N — 5 the emission 
intensity is. 



200e"^^r - BOe^^^r. (C5f) which is shown in Fig. [4(a)] 



I^{t) = -[162e"^^-|-16e"^^(24T-l) + e~^^(240T--143)], 

(C6) 
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